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The properties of color-flavor locked strange quark matter in an external strong magnetic field are investigated 
in a quark model with density-dependent quark masses. Parameters are determined by stability arguments. It is 
found that the minimum energy per baryon of the color-flavor locked (MCFL) matter decreases with increasing 
magnetic-field strength in a certain range, which makes MCFL matter more stable than other phases within a 
proper magnitude of the external magnetic field. However, if the energy of the field itself is added, the total 


energy per baryon will increase. 
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I. INTRODUCTION 


Strange quark matter (SQM) is an interesting topic [1] not 
only because of its great theoretical significance, but also its 
many applications, e.g., in studying quantum chromodynam- 
ics (QCD) phase diagrams [2-6], properties of strangelets [7, 
8], and the structure of compact stars [9, 10], etc. In 1984, 
Witten conjectured that quark matter with strangeness might 
be the true ground state of QCD [11]. Soon after Witten’s 
conjecture, Farhi and Jaffe studied the stability of SQM with 
the conventional MIT bag model and found that SQM is ab- 
solutely stable around the normal nuclear density for a wide 
range of model parameters [12]. Since then, SQM has be- 
come a main topic in a number of meaningful works [13—20]. 


It has been demonstrated that SQM at high density may be 
in the color-flavor locked (CFL) phase where quarks with dif- 
ferent color and flavor quantum numbers form Cooper pairs 
with a large binding energy [21]. It is, thus, possible that CFL 
matter, rather than nuclear matter, may be the ground state of 
strange quark matter at high density. Therefore, a compact 
star is suggested to include color superconducting quark mat- 
ter in its inner core [22]. 


It is generally believed that properties of quark matter have 
been strongly affected in the presence of a strong magnetic 
field [23, 24]. A strong magnetic field widely exists on 
the surface of stars. The observed magnetic field strength 
on the surface of pulsars could be in the order of 10!?— 
10'°G. And the magnetic-field strength on the so-called 
magnetars could be in the order of 10'4-10!°G or even 
higher [25, 26]. In fact, the biggest magnetic field that can 
be sustained by strange stars was estimated to be as large as 
1.5 x 107° G [27]. Although the origin of the strong magnetic 
field is not completely clear and still under active investiga- 
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tions, some ways to understand its existence are avaible, e.g., 
the amplification of the relatively small magnetic field during 
the star’s collapse with magnetic flux conservation [28],and 
the magneto-hydrodynamic dynamo mechanism with large 
magnetic fields generated by rotating the plasma of a pro- 
toneutron star [29]. Furthermore, noncentral high-energy 
heavy-ion collisions could generate intense magnetic fields 
as high as about 10/9 G [27], corresponding to eBm ~ 6m?2, 
where e is the fundamental electric charge and m, is the pion 
mass. It is therefore necessary to study the properties of CFL 
in the presence of an external magnetic field. 


In past years, magnetized strange quark matter (MSQM) 
and CFL matter have been studied with many phenomeno- 
logical models, e.g., the bag model [30], the Nambu-Jona- 
Lasinio (NJL) model [31-36], and the mass-density depen- 
dent model [37, 38], etc. The MCFL matter has a wide 
range of model parameters characterized by the so-called sta- 
bility window [39], and has also been studied in the NJL 
model [40—44], as well as in the quasiparticle model [45, 46]. 


As is well known, particle masses vary with environment, 
i.e., they depend on density or chemical potential. The equiv- 
particle model [1, 2] takes this effect into account by density- 
dependent quark masses. In recent years, this model has been 
extensively applied to study the properties of SQM [1, 15, 47— 
55]. In this paper, we extend it to investigate the properties of 
CFL matter when a strong magnetic field appears. It is found 
that MCFL matter is more stable than the other phases within 
a proper magnitude of the magnetic field. At a fixed density, 
the energy density of MCFL matter varies with the magnetic 
field strength. At Bm > 10! G, the energy per baryon, pres- 
sure, and quark chemical potentials get smaller because the 
quantum number of the corresponding Fermi momentum ap- 
proaches to zero. 


This paper is organized as follows. In Section II, we give 
the thermodynamic formulas used for the study of MCFL 
matter in the equiv-particle model with density-dependent 
quark masses. Then we present the numerical results and dis- 
cussions in Section III. Section IV is a short summary. 
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Il. THERMODYNAMIC TREATMENT 


Our starting point is the thermodynamic potential density 
of a free-particle system, i.e. 


29: 
apip (F) o 


where the summation index, i, goes over u, d, s quarks and 
electrons, m; is the corresponding particle mass, ju; is the 
chemical potential, g; is the degeneracy factor with a value of 
3 for quarks and 1 for electrons, while the degeneracy due to 
spin has been denoted by a factor of 2. 

In the case of CFL phase, due to the energy gap, A, de- 
termined by solving the gap equation, a new term should be 
added to the above expression. The thermodynamic potential 
density of CFL matter is then 


29; 3A? 7? 
fom = GBs | (&- m) Bp o +B, ®© 


where €; = \/p? + m? is the dispersion relation of a free par- 
ticle with a mass of mj, i = (Hu + Ha + Us) /3 is the average 
of the quark chemical potentials. The second term is from the 
pairing contribution, and the last term, B, is the famous MIT 
bag constant to take the vacuum energy into account. 

To consider the effect of a magnetic field, we assume a 
constant magnetic field with a strength of Bm along the z 
axis. Due to Landau diamagnetism, the single particle energy 
spectrum can be written as 


Eil = V p22 +m? + 2\qi| Bm[l + 1/2 —sen(qi)S], (3) 


where p, is a component of the particle momentum along the 
direction of the magnetic field, q; is the electric charge of 
quarks when qu = 1/3, qa = qs = —1/3, qe = —1, l = 
0,1,2,... is the principal quantum number for the allowed 
Landau levels, and S = +1/2 refers to spin-up and spin- 
down states, respectively. The sign function ‘sgn’ equals 1 
with a positive argument and —1 with a negative argument. 
For the sake of convenience, one normally sets v = l + 
1/2 — sgn(q;)S, where v = 0,1,2,.... The single particle 


energy becomes 
Eivy = [P+ M}, (4) 
where M; y = ym? + 2v |qi| Bm. 


The integration over the p,—p, plane in the momentum 
space should be replaced by 


J | dpsðpy > 2r lail Ba S902 0). (5) 


After this substitution, Eq. (2) becomes 


Pern =D | (Ei, — i) dpz 


2a 2 


+ B, (6) 
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where we have used the notations: 


gi |Gi| Bm 
fay = ltl Bagg 


Piv =f wi? — M2, a 


— Jv,0): (7) 


The upper bound Vmax of the summation index v in Eq. (6) 


is 
2 pe 
Vmax = int Gas (8) 
2 gil By 


where the function int(x) means taking the integer part of its 
argument z. 

After carrying out the integration in Eq. (6), we have 
the thermodynamic potential density in the conventional bag 
model 


Vmax 
NRucrr =- > Y fiw e pi? — MZ, 
i v=0 
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+ B. (9) 


To include medium effect, the quark masses should be den- 
sity/chemical potential dependent. In the chemical potential 
dependent case, one can use the quasiparticle model, as has 
been done in Ref. [45]. In the density-dependent case, the 
actual chemical potential, j:;, should be replaced with an ef- 
fective chemical potential, yw; [1], i.e. 


34202 
m= a- — +B, (10) 


where j is now understood as the average of the effective 
chemical potentials, and 2; is connected to the effective 


chemical potentials by 
Vmax už 
Mix) | 


a 5 fiw È 4/ we? — MŽ, = M? „arcsh ( 
v=0 
(11) 


All other thermodynamic quantities can be derived from 
No. Specially, the particle number density is given by n; = 
—0No/Om;, giving 


Vmax 


i|di| Bm 2A? i 
iG a >> (2-b0)4/ uf? — M2, +=". (12) 
v=0 


The energy density for the MCFL matter is then 


Emori = 20+ 5 [Lj Ni. (13) 
Upon application of Eqs. (10)-(12), we have 


Emcri = > sit a Ee ofui yu — MŽ, 


i 3A2 72 
os Hi )| wie , +B. (14) 
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Because the quark masses are density dependent, the actual 
chemical potential is generally not equal to its effective one. 
In fact, they are linked by 

Hi = My + pi, (15) 
where uy is due to the density dependence of the quark 
masses. Its explicit expression is obtained by applying the 


fundamental differential equality dEMCFL = ye idn; with 
Eqs. (10)-(15), giving 


ON dmj 
p= 2 Om; dn; , 


9;|¢;|BmD 
7 1842n4/3 (2 — ôvo)mjarcsh 
Uv 


(16) 


Lj 
(17 
E ) 


Because electrons do not participate in strong interactions, 
their actual chemical potential is equal to the effective one, 
i.e., le = ple. 

Due to the external magnetic field, the longitudinal pres- 
sure and transverse pressure become different, i.e. 


Py = -9o + pani (18) 


and 


L = — 2+ $ uni — MiBm, (19) 


where P) is the total parallel pressure and P, is the transverse 
pressure. The system magnetization is given by 


AN Qi 82; dMiy 
ie OBn D (+ + Mi, dBm ) s 


i 


Upon application of Eqs. (10), (11), and (20), we have the 
following explicite expressions 


= a pon ae 


2 D 
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«2 2 * 
fe My +i 3420? 
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HI. THE PROPERTIES OF MCFL MATTER 


In the equiv-particle model, the quark mass can be devided 
into two parts as 
Mi = Mio + mM, (23) 
where mio is the quark’s current mass, and my represents the 
effect due to the interaction between quarks. In principle, the 
density dependence of my should be determined from QCD. 
As mentioned before, however, there is no way to exactly 
solve QCD presently. Therefore, the density dependence is 
normally given phenomenologically. It can be shown that the 
following parametrization is reasonable, 
(24) 


Mi = Mio + Pal 1/3 ’ 
where D is a fixed constant determined by stability argument, 
n is the total baryon number density, and the exponent of the 
baryon number density was derived based on the in-medium 
chiral condensates and liner confinement at zero temperature. 
Such a form satisfies lim,_,97m, = oo and limps mr = 0, 
which are the requirements of quark confinement and asymp- 
totic freedom. 

Because weak equilibrium is always reached in SQM, rel- 
evant chemical potentials satisfy 


Md = Hs, (25) 


and 


Hu + He = Hs- (26) 


Therefore, the effective chemical potentials also meet the 
correspongding relations 


ba, + He = Ma (27) 
and 
ee. (28) 
We also have the baryon number density 
1 

n= g (mu + na + ns), (29) 

and the charge density 

1 1 

Q = gnu 3 itd gns ne. (30) 


The charge-neutrality condition requires Q = 0. 

For a given total baryon number density of n, we can obtain 
the respective u}, už, HŽ, and ue by solving Eqs. (27)-(30) 
with the help of Eq. (12). The number densities ny, ng, ns, 
and ne can then be obtained. The energy density is calculated 
by Eq. (14), while the pressures are, respectively, obtained 
from Eqs. (21) and (22) for different values of B, D, and By. 

In Fig. 1, we give the energy per baryon of SQM, MSQM, 
CFL matter, and MCFL matter, respectively, as a function of 
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Fig. 1. Comparison of the energy per baryon versus the baryon num- 
ber density for SQM, MSQM, CFL matter, and MCFL matter, re- 
spectively. The pressure at the energy minimum is exactly zero for 
each quark matter phase. 


the baryon number density. Because the current mass of u/d 
quarks are very small, we simply take muo = Mao = 0. But 
for the current mass of the strange quark, we take mso = 
80 MeV. The electron does not participate in the strong in- 
teraction, its mass is very tiny, and also ignored. For conve- 
nience of comparison with previous works, we take the paring 
parameter to be A = 100 MeV. The bag constant, B, and the 
confinement parameter, D, are taken to be B!/4 = 140 MeV 
and D!/2 = 120MeV. In our calculation, we assume the 
magnetic field to be constant with its direction along the z 
axis. Because the system will become unstable when the 
magnetic field strength is higher than 10?° G, as discussed 
by Chakrabarty [23], we take the magnetic field strength to 
be Bm = 10/°G. From top to bottom in Fig. 1, there are 
three features. Firstly, the energy minimum (the solid trian- 
gle) corresponds exactly to the zero pressure (open circle) for 
each case. In fact, the exact coincidence of the lowest energy 
state and zero pressure is a basic requirement of the funda- 
mental thermodynamics, as pointed out in Ref. [15], and de- 
rived in detail in Ref. [1]. Secondly, the energy per baryon of 
CFL matter and MCFL matter is lower than that of SQM and 
MSQM. So we can see that the quark pairing effect greatly 
increases the stability chances of SQM. Thirdly, the energy 
per baryon of MSQM and MCFL matter is lower than that of 
SQM and CFL matter, respectively. We can see that the ex- 
ternal magnetic field in a proper magnitude lowers the energy 
per baryon through the rearrangement of the Landau energy 
level of magnetized quark matter. Generally, we have the in- 
equality relation of the energy per baryon as 


E E E E 
= —— < — — 


"imort “@lcrnr ™IlmMsqm ” 


(31) 


SQM 


In Fig. 2, we give the minimum energy per baryon of 
MCFL phase as a function of the magnetic field strength. 
When the magnetic field strength, Bm, is small, it is obvi- 
ous that the energy is nearly constant. The energy per baryon 
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Fig. 2. The minimum energy per baryon of MCFL matter as a func- 
tion of the magnetic field strength. Relevant parameters are indicated 
in the figure. 


starts to decrease obviously as a function of the magnetic 
field between 1018 G and 10!° G. When the magnetic field 
strength exceeds 101° G, the energy per baryon decreases 
quickly. Therefore, an external magnetic field with proper 
strength lowers the energy per baryon. In this regard, one 
should note that the energy from the external magnetic field 
was not added. Otherwise, the total energy per baryon will 
increase. 

There are different views on whether or not the energy con- 
tribution from the magnetic field should be included. If one 
would like to include the field contribution, one should know 
how the quark matter produce the magnetic field. As men- 
tioned in the introduction, the origin of the the strong mag- 
netic field is presently not very clear, although some ways to 
understand it are available. Therefore, we treat the magnetic 
field as an externally forced field. 


P (MeVifm*) 


m,,= 80 MeV 


D"?=160 MeV 
AA =100 MeV 


10° 10” 10° 10° 10” 


Fig. 3. The pressure in MCFL matter as a function of the magnetic 
field strength at two (the solid line) and three (the short-dot line) 
times the nuclear saturation density no = 0.165 fm™?. 
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As is well known, the space becomes anistropic when an 
external magnetic field is presented. To compare the mag- 
nitude of the longitudinal and transverse pressures, we plot, 
in Fig. 3, the Pj and P, at the given densities of n = 2n 
(the solid line) and n = 3npo (the short-dot line) where 
no = 0.165 fm~? is the nuclear saturation density, as func- 
tions of the magnetic field strength. The difference between 
P| and P, reflects the breaking of the rotational symmetry by 
the magnetic field. We can see that the pressure stays as con- 
stant when the magnetic field strength is lower than 10!°G. 
When the magnetic field strength is larger than 10'°G, the 
pressure anisotropy starts to become noticeable: the parallel 
pressure, Pi, increases far beyond the constant value, while 
the vertical pressure, P, , decrease from the constant value. 
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Fig. 4. The chemical potential of quarks in MCFL matter as a func- 
tion of the magnetic field strength for n = 3no. 


In Fig. 4, the chemical potentials are shown as functions of 
the magnetic field strength for n = 3no. When the magnetic 
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